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Abstract 

Properties of cold nuclear matter are studied within a generalized Nambu-Jona-Lasinio model 
formulated on the level of constituent nucleons. The model parameters are chosen to reproduce 
simultaneously the observed nucleon and pion masses in vacuum as well as saturation properties 
of nuclear matter. The strongest constraints on these parameters are given by the empirical values 
of the nucleon effective mass and compression modulus at nuclear saturation density. A preferable 
value of the cut-off momentum, determining density of active quasinucleon states in the Dirac 
sea, is estimated to about 400 MeV/c. With the most reasonable choice of model parameters we 
have found a first order phase transition of the liquid-gas type at subsaturation densities and the 
gradual restoration of chiral symmetry at about 3 times the saturation density. Fluctuations of the 
scalar condensate around its mean-field value are estimated and shown to be large in the vicinity 
of chiral transition. 

PACS numbers: ll.30.Rd, 12.39.Fe, 21.65. +f 
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I. INTRODUCTION 



Properties of strongly interacting matter, especially its possible phase transitions, are 
presently in the focus of experimental and theoretical investigations. At zero baryon density 
and high temperature the situation becomes more and more clear due to continuous progress 
in lattice calculations. Latest simulations jl| including dynamical quarks with realistic 
masses predict a crossover-like deconfinement (chiral) transition at temperatures around 
170 MeV. The situation at finite baryon densities remains uncertain and subject to model 
building. Different models predict many interesting phenomena in this case, ranging from 
restoration of chiral symmetry to color superconductivity. 

In our recent analysis ^ we tried to find phenomenological constraints on possible phase 
transitions at high baryon densities. The main constraint comes from the very existence 
of a bound state of symmetric nuclear matter at low temperature and baryon density 
p ^ 0.17 fm~^ . We have found that predictions of some popular models are in direct con- 
tradiction with this fact. For example, several quark models predict new phases which have 
too low pressure as compared to ordinary nuclear matter with the same baryon chemical po- 
tential. Therefore, such phases are unstable with respect to hadronization. Our aim in this 
paper is to construct a QCD motivated effective model which is able to reproduce correctly 
the nuclear saturation point and then examine its predictions at higher densities. 

Many phenomenological models have been suggested to describe properties of cold nuclear 
matter in terms of nucleonic degrees of freedom. The relativistic mean-field models of the 
Waleoka type S Q tun.d out to be .ath. successful iu dcscnbiug p.opct.cs of .cduun aud 
heavy nuclei. One should mention also microscopic models, dealing with nucleon-nucleon 
forces. Parameters of these forces are usually extracted from scattering data and the observed 
binding energies of light nuclei (see e.g. Ref. P| and references therein). In addition to their 
numerical complexity, such approaches are essentially nonrelativistic and therefore, their 
accuracy diminishes at high densities. All these models have a serious drawback, namely, 
they do not respect the chiral symmetry of strong interactions. 

There exist several chiral models which could potentially be used for describing nuclear 
matter. Most popular are the Nambu-Jona-Lasinio (NJL) model and the linear sigma 
model ~\. They are able to explain spontaneous breaking of chiral symmetry in vacuum 
and its restoration at high energy densities. But the simplest versions of these models fail to 
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reproduce nuclear saturation properties. In particular, the linear sigma model predicts only 
abnormal state of nuclear matter jjsl where the chiral symmetry is restored and nucleons 
have vanishing effective mass. Several more complicated models of this kind have been 
suggested in Refs. 0, 10, 12, [l^. Although they are able to reproduce the nuclear ground 



state, some new problems appear within these models. In particular, some of them do not 
predict restoration of chiral symmetry at high baryon densities. 



ter 



There were also attempts to use the NJL model for describing cold nuclear mat- 
] . It was argued in Refs. ^ , Q| that bound nucleonic matter with spon- 
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taneously broken chiral symmetry is not possible within the standard NJL. This conclusion 
was based on implicit assumption that the maximum (cut-off) momentum for constituent 
nucleons A is the same as normally used in quark models (A ~ 0.6 GeV). The authors 
of Ref. [l^ suggested to include additional (scalar-vector) interaction terms to reproduce 
Dserved saturation properties of nuclear matter. On the other hand, it was shown in 



the o 

Ref. [1^ that by assuming sufficiently low value of cutoff-momentum (A ~ 0.3 GeV) it is 
possible to achieve a bound state at normal density even in the standard NJL model. How- 
ever, in this case the nucleon effective mass at p = po is predicted by a factor of two smaller 
as compared to its empirical value. 

Below we reconsider the possibility of using the NJL model for describing cold nuclear 
matter. We use the generalized version of this model including additionally the scalar-vector 
interaction. We also take into account explicit symmetry breaking effects by introducing 
a bare nucleon mass. Indeed, with this modifications a good agreement with observed 
saturation properties of nuclear matter is achieved. 

The paper is organized as follows. In Sect. II we formulate the model and explain the 
procedure of choosing parameters. Comparison of our results with observed characteristics 
of nuclear matter is given in Sect. III. In Sect. IV we analyze characteristics of the chiral 
symmetry restoration and demonstrate the importance of chiral fluctuations. In Sect. V we 
summarize and discuss the results. 

II. FORMULATION OF THE MODEL 

I— I 

In this paper we return to the original formulation of the NJL model jS] in terms of 
nucleonic degrees of freedom. Namely, we consider a system of nucleons, interacting via 
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point-like scalar (pseudoscalar) and vector (axial vector) 4-fermion terms. The correspond- 
ing Lagrangian for a spinor nucleon field ip is written as {h = c = 1) 



where mo is the "bare" nucleon mass, r are isospin Pauli matrices, are standard Dirac 
matrices and 75 = i7°7^7^7'^ . Gs and Gy are the coupling constants for the scalar (pseu- 
doscalar) and vector (axial vector) NN interactions, respectively. 

As will be clear later, this Lagrangian is too restrictive to reproduce saturation proper- 
ties of nucleonic matter at zero temperature. In particular, the model should explain the 
phenomenological fact that binding energy per nucleon of homogeneous isospin-symmetric 
matter is bound by about 16 MeV at the equilibrium density p = po — 0.17 fm~'^ . In other 
words, this means that at this density the energy per baryon reaches its minimum value 

E 
B 



TUN - 16 MeV , (2) 

P=PO 



where ttt-at = 938 MeV is the nucleon mass in vacuum. 
Following Refs. 



5. below we consider a generalized version of the NJL model including 

an additional 8-fermion scalar-vector (SV) interaction term. The full model Lagrangian is 
written as 



C = Co + ^ - (^75T^)'] [(^7^^)' + (^7m75^)'] . (3) 

It is easy to see that at mo = this Lagrangian is chirally symmetric, i.e. it is invariant with 
respect to SU{2)l ® SU{2)ji transformations. Below we investigate sensitivity of nuclear 
matter equation of state to the coupling constants Gs, Gy and Gsv ■ In fact, they will be 
chosen to fulfil the condition (0). 

Within the mean-field (Hartree) approximation we replace {ipTip)'^ operators by 
2il)Til) < ipTip > — < ipTip >^ , where F is any combination of matrices appearing in the 
interaction terms of the Lagrangian (jH)). The angular brackets denote quantum-statistical 
averaging at fixed temperature T and baryon density p . In the considered case of isospin- 
symmetric matter at T = all terms containing 75 vanish and the Lagrangian takes the 
form 



4 



'Cmfa = ipiic^-m- 7°Sv/) ip -U{p, ps) 



(4) 



where m, Sy and U are functions of scalar {ps =< "ipip >) and vector (p =< ijj'y'^ip >) 
nucleon densities. These quantities are defined as follows 



m = mo - Gsps , 

Sy = (Gy - GsvpDp, 

U{p, Ps) = \ {GspI - Gyp^ + SG^ypVI) • 



(5) 
(6) 
(7) 



In Eq. (0) we have introduced a renormalized scalar coupling 



Gs = Gs^ Gsvp" = Gs[l + a [p/ p,f] , 



(8) 



where a is a dimensionless constant, a = GsvPo/Gs ■ Solutions of the gap equation 
denoted below have the meaning of the nucleon effective mass. 

The Euler-Lagrange equation for the fermion field ip has the form of Dirac equation with 
constant effective mass m and vector potential Sy. By using the plane wave decomposition 
of ip one gets the explicit relations 



where i^a? = 4 is the spin-isospin degeneracy factor for nucleons, nN{p) and nj^{p) are the 
Fermi-Dirac occupation numbers for nucleons and antinucleons. The last term (—1) in 
square brackets of Eq. ^ originates from the negative energy levels of the Dirac sea. It 
gives divergent contribution which needs to be regularized. Following common practice we 
do this by introducing 3-dimensional cut-off momentum A , i.e. replacing —1 by — 0(A — p) 
where G(x) = (1 + sgnx)/2 . One can interpret A as maximum momentum determining the 
number of active quasinucleon levels in the Dirac sea |2^. 




From Lagrangian one can easily derive the energy-momentum tensor. The energy 
density of homogeneous nuclear matter is written as 




(10) 





e = 




- + rn? [uMij)) + n^(p) - 0(A - p)] + Syp + U (p, ps) + eo 
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where again the three-dimensional regularization is apphed to the divergent momentum 
integral. The constant eo is introduced to set the energy density of the vacuum state (p = 0, 
m = m^v) equal to zero. Obviously, energy per baryon is E/B = tj p . 

Below we consider isospin-symmetric matter at zero temperature. In this case njj^^p) = 
and riNip) = 0(pf — p) where pf is the Fermi momentum related to the baryon density, 



6n^ 



(12) 



The scalar density reads 



A 



PS = -J^N 

PF 



d^p m uj^m 



(27r)3 ^p2 + ^2 47r2 



(13) 



where 



$(x) = yrT^-^in ^^zji^^; . (14) 

^ V 1 + — 1 

At sufficiently low baryon densities corresponding to pp < A the integral in the r.h.s 
of Eq. ()13|1 is positive and ps < ■ However, at high densities, when pp > A, the in- 
tegral becomes negative and, therefore, > . Then, according to Eq. ^ the nucleon 
effective mass becomes smaller than mo . Of course, at large enough momenta the approx- 
imation of a point-like NN interaction will fail due to the nontrivial internal structure of 
nucleons. In this situation one should explicitly consider quark degrees of freedom. 
Expressing ps iii terms of m from Eq. (0), we rewrite Eq. as 

A 

VP + H \ h eo . (15) 



2 , Vl +x2 + 1 



PF 

One can consider this expression as a functional of two independent variables, m and p . Then 
the gap equation follows from the minimization of e (m, p) with respect to m . As has been 
already noticed in Ref. Q|, the inclusion of SV interaction is equivalent to renormalizing 
the scalar coupling, Gs — ^ Gs ■ 

From thermodynamic identities at T = 0, 

one can calculate the baryon chemical potential p and pressure P as functions of baryon 
density p. It is obvious that pressure should vanish at the saturation point, i.e. -P(po) = . 
This gives a nontrivial constraint on the model parameters. 



Our model contains five adjustable parameters, A, tuq, Gs, Gy, a, which we fix by fit- 
ting the vacuum masses of nucleons and pions, as well as the saturation properties of 
nuclear matter. First constraint follows from the gap equation in the vacuum. Taking 
p = 0, Pf = 0, m = mjq in Eqs. (0)) (fT^ one has 

= mo - GspT = + ^nGs"^^ (^) . (17) 



Another constraint [16j comes from the soft pion phenomenology. By using the La- 
grangian (Q) in the random phase approximation one can calculate the polarization operator 



for pionic excitations (for details see Ref. 19[). By applying further the PCAC equation for 



the axial-vector current of nucleons in vacuum one gets 

mlfl = m,\pT\. (18) 

where ~ 140 MeV and ~ 93 MeV are the vacuum values of pion mass and pion decay 
constant, respectively. This equation is valid in the lowest order approximation in m,r/A 
and can be considered as analogue of the Gell-Mann-Oakes-Renner relation [2^ following 
from the quark structure of mesons. According to Eq. ()18p the pure chiral limit mo = , 
used e.g. in Ref. Q|, formally corresponds to vanishing pion mass. 

Choosing different cut-off momenta A as input we use Eqs. (jl7p -(|18| ) to find mo and Gs- 
The parameters Gy and a can be fixed now by the requirement that e/p attains the minimum 
value (0) at p = po • At fixed A we first find m = m*j^{p) from the gap equation (see Eqs. (jSJ, 
()13|)) and then calculate the energy per baryon by using Eq. (fT3j) . The calculation shows 
that at A ^ 0.2 GeV it is not possible to reproduce the saturation point with observable 
binding energy and density at any a . This becomes possible only at A ^ 0.3 GeV 

An illustration of the procedure used to determine parameters Gy and a is given in Fig. 
for the case A = 0.4 GeV. At fixed a we choose the ratio = Gy/Gs to obtain minimal 
energy per baryon given in Eq. Q. As one can see in this figure, increasing a shifts the 
minimum in e/p to smaller baryon densities. In the considered case the choice a = leads 
to abnormal (m^ ~ 0) bound state at density p ~ 2.2po . The correct position of the 
saturation point is obtained for a = 0.032 . We repeat this analysis for other values of A 
from the interval 0.3 GeV < A < 0.7 GeV . 

The corresponding parameters are listed in TableHl The second and third columns give Gs 
and mo values, determined, respectively, from Eqs. (fTTj) and fTH|). 
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FIG. 1: Energy per baryon in cold nuclear matter as a function of baryon density p, calculated 
with the cut-off momentum A = 0.4 GeV. Different curves corresponds to different values of scalar- 
vector coupling a . The parameter ^ denotes the ratio of vector and scalar coupling constants. 



TABLE I: Parameters of the NJL model 



A (GeV) 



Gs(GeVfm^ 



mo (MeV) 



a 



0.3 
0.4 
0.5 
0.6 
0.7 



3.653 
1.677 
0.900 
0.542 
0.354 



95.7 
41.3 
21.7 
12.9 
8.41 





0.032 
0.026 
0.024 
0.0225 



0.885 
0.931 
1.226 
2.224 
4.506 



The last two columns show the parameters a and ^ which give the best fit of the saturation 
point at fixed cut-ofT momentum A . One can see from the table that the parameter a 
slowly decreases with A at A ^ 0.4 GeV. On the other hand, the relative strength of vector 
interaction ^ grows noticeably with increasing A . It is natural to assume that the bare 
nucleon mass should not be smaller than Sniog , where mog = (5±1) MeV 2l| is the isosp in- 
averaged current mass of light quarks. From this point of view the parameter sets with 
A > 0.6 GeV should be excluded. Due to the same reason the parameter set with mo = 
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and 
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FIG. 2: Energy per baryon as a function of p for different A. 



A = 0.641GeV, G5 = 0.459 GeVfm~^ a = 0.023 , ^ = 2.92 
suggested in Ref. [l^, does not seem reasonable. 



(19) 



III. NUCLEAR MATTER WITHIN THE NJL MODEL 



Now we present the predictions of the generahzed NJL model formulated in the preceding 
section. Figure El shows density dependence of energy per baryon calculated with parameter 
sets from Table HI Although the results are similar at p ^ po ? at higher densities one can 
see strong sensitivity to A. The curve for A = 0.3 GeV shows especially large deviation from 
the others. As already noted above, one should be very cautious by applying the NJL model 
at high densities, corresponding to pp > A . This is particularly important for sets with 
smaller A . Indeed, for A = 0.3 GeV pp > A holds already at p > 1.4po 3^ . 

Figure El shows density dependence of the nucleon effective mass m*^ calculated with the 
same parameter sets. On can see that the behaviour depends strongly on A even at low 
densities. For example, at A = 0.7 GeV it has a very peculiar behavior. In this case 
first increases, reaching maximum value m*^ ~ 1.4 GeV at p ~ lOpo , and then drops down 
approaching the bare mass at p ~ 17po. For all A the nucleon effective mass becomes small 
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FIG. 3: Nucleoli effective mass as a function of baryon density for different A. The dots show 
estimates based on the QCD sum rules. 

at sufficiently high densities, but it does not vanish. This would be possible only in the case 
of exact chiral symmetry (mo = 0) when = is an exact solution of the gap equation. 
Therefore, for all parameter sets with tuq ^ the chiral symmetry is only approximately 
restored at high densities. By filled dots in Fig.Elwe show the mass values obtained from the 
QCD sum rules 3]- Although the error bars are large, these estimates clearly favour the 
parameter set with A = 0.4 GeV. As compared to vacuum, a 30% reduction of the nucleon 
mass at p = po is predicted in this case. The results for different A are presented in Table ITTl 

TABLE II: Properties of nuclear matter within the NJL model 



A (GeV) 



mlf{po)/mN 



Kipo){MeY) 



S^^(MeV) 



'0 



0.3 
0.4 
0.5 
0.6 
0.7 



0.385 
0.683 
0.834 
0.905 
0.944 



1350 
285 
322 
347 
337 



101 
42.1 
25.1 
15.9 
10.9 



1.24 
2.66 
5.08 
9.51 
16.3 



The compression modulus of nuclear matter, K, is an important characteristics of the 
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FIG. 4: Compression modulus as a function of baryon density at different A. 



equation of state. It is defined as 



K dP 



9 dp ^ dp ^ 



d'e 



dp 



2 ' 



(20) 



where the second and third equalities follow from Eq. (fTH|l . The density dependence of K 
is shown in Fig. |3] for different values of A . Negative K values correspond to mechanical 
instability of matter with respect to density fluctuations. The density interval with i^' < 
is known as the spinodal region of a first order (liquid-gas) phase transition. As one can see 
from Fig. ^ this phase transition is undoubtedly predicted in all cases considered. This is 
not surprising because such a phase transition is a consequence of the saturating property 
of the equation of state. Indeed, since pressure is zero at p = and p = po is must be 
a nonmonotonous function of p . Fig. 0] also shows that the compression modulus changes 
very rapidly in the vicinity of the saturation point, especially in the case of A = 0.3 GeV. 
The values of i^' at p = po are given in Table |nj One can see that the parameter set 
with A = 0.4 GeV give s the K value which is close to the phenomenological estimate, 
is: ^ 250 - 270 MeV [2^. 

In studies of in-medium effects, the so-called pion-nucleon sigma term, Sttat , is often 
used • At small densities it can be defined as 

s.iv = ( 1 — : r : i ^ , (21) 



< qq >vac 



p 
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FIG. 5: Pion-nucleon sigma term as a function of baryon density at different A. 
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FIG. 6: Energy per baryon calculated within different nuclear models (see text). 



where < gg > is the quark condensate or, in the context of our model the quark scalar 
density. Assuming that the scalar nucleon and quark densities are proportional to each other, 
we calculate Sttat by using Eq. (j7T|) with the replacement <qq> / <qq >vac— ^ Ps/Ps"^- The 
results of this calculation are shown in Fig. El The S^^r values at p — > are also given in 
Table ITTl At A = 0.4 GeV the calculated value is close to the estimate Sttat — (35 ± 5) MeV, 
obtained from the chiral perturbation theory 



12 



Based on the results presented above we come to the conclusion that the generalized NJL 
model with A ~ 0.4 GeV gives the best description of cold nuclear matter in comparison 
with other choices of A . In Fig. [Ul we compare energies per baryon predicted by our model 
with the results of other nuclear models namely, the Walecka-type TMl model 

Q, the 

Chiral Hadron Model (CHM) jisl and the Variational Chain Model (VCM) Note that 
the VCM includes two- and three-body NN forces and takes into account the correlation 
effects. One can see that at densities p ^ 2 po out model and other models give very similar 
results. Thus, with A = 0.4 GeV our fit of nuclear saturation properties is at least of the 
same quality as in other models listed above. 



IV. IDENTIFICATION OF THE CHIRAL TRANSITION 

In this section we examine properties and possible signatures of the chiral transition 
in cold nuclear matter. As noted in Ref. j23] the direct signature of the chiral symmetry 
restoration is not the effective mass (the latter may even increase with energy density), but 
the vanishing chiral condensate. In our approach the chiral condensate coincides with the 
scalar nucleon density ps =< ipip > ■ Figure [7| displays this quantity calculated for the same 
parameter sets as before. Unlike the effective mass, \ps\ decreases practically linearly with p 
at not too high densities. We define the critical baryon density of the chiral transition, pc, 
by the linear extrapolation of ps{p) to zero. Numerical values of these densities (in units 
of po) are given in the last column of Table ITTl 

The results of Fig.[7|suggest that scalar density curves for different A have similar behavior 
and roughly can be obtained from one another by rescaling the axes. In Fig. |Hlwe show the 
same results, but scalar and vector nucleon densities are now scaled by the factor A~'^ . 
One can see that in these variables the results for different A are indeed very close to each 
other. It is evident now that the critical density for the chiral transition may be roughly 
estimated as 

Pc ~ 0.06A3 . (22) 

The corresponding Fermi momentum p^r at p ~ Pc equals approximately 0.96A. Therefore, 
the region of the chiral transition is close to the boundary of the model applicability (see 
discussion above). 

It is instructive to compare the vacuum value of nucleon scalar density with the quark 
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FIG. 7: Scalar condensate density as a function of baryon density for different A. 
8 
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FIG. 8: Same as Fig. [71 but with scaled densities, 
condensate <qq>. The phenomenological value of the quark condensate in vacuum is 



I < gg > l^ac ~ 2 I < MM > l^ac ^ 2 (225 ± 25 MeV)^ ^ (3 ± 1) fm"^ . (23) 

According to Fig. [71 for A = 0.4 GeV one gets the value \ps^^\ — 0.54 fm~^ which is by about 
factor 6 smaller than \ <qq > |vac- For A = 0.5 GeV the ratio is close to 1/3 which might 
be expected from the naive constituent quark model. For higher A the predicted nucleon 
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FIG. 9: Pressure as a function of chemical potential for different A. The curve labelled by VCM 
is calculated within the model of Ref. j^. The insert shows -P(//) in the region of liquid~gas phase 
transition. 

condensate is comparable with < qq >vac • 

When looking for possible phase transitions it is most useful to calculate pressure P as a 
function of chemical potential /i . A first order phase transition is signalled by appearance 
of several branches of P{fi) ■ At fixed fi only the branch with highest pressure corresponds 
to a stable phase. The point of phase transition, /x = /Xc , is given by intersection of a 
stable and a metastable branches. At this point dP/dfi exhibits a jump which corresponds 
to difference of densities between coexisting phases. Figure E] shows the pressure curves 
predicted by our model for different A . At small fi one can see the presence of a first-order 
phase transition which takes place at fi c:^ rriN — 16MeV ^ 922 MeV for all considered sets 
of model parameters. The region of mixed phase covers subnuclear densities, p < po . This 
is a liquid-gas phase transition which is well-known in nuclear physics. It is responsible for 
the multifragmentation phenomenon in hot nuclear systems. There is no clear evidence of 
other non-trivial behaviors of P{p) ■ Calculating P as a function of p shows that pressure 
varies rather smoothly at p ~ pc i-e. in the vicinity of the chiral transition. 

The results presented above were obtained assuming validity of the men-field approxi- 
mation. In particular, fiuctuations of the scalar condensate ps or, equivalently, the nucleon 
effective mass m around their values determined from the gap equation were completely 
disregarded. As we shall see below, this approximation becomes rather poor at p ~ • To 
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FIG. 10: Energy per baryon as a function of nucleon mass at fixed baryon densities (shown in units 
of po in the key box) within the NJL model with A = 0.4 GeV. 

estimate possible magnitude of mass fluctuations, it is instructive to study the energy per 
baryon, e/p, as a functional of p and m, i.e. without using the gap equation. By fixing 
baryon density, one obtains the profiles of e/p as functions of m . Their minima correspond 
to m = m*j^{p) . The results of calculations are shown in Fig. El for the parameter set with 
A = 0.4 GeV. 

As seen in Fig. at p — energy per baryon is a sharp function of mass near its 
minimum at m ^ mTv. As p grows, the minimum of e/p shifts to smaller masses and 
its width significantly increases. At p = po the minimum corresponds to the saturation 
point where m = ~ 0.683 mAr. When p approaches the characteristic density of chiral 
transition p ~ p^ the minimum shifts to m ~ and e/p curves become rather fiat. For 
example, at p = 2po the width of mass distribution becomes comparable with the value 
m*^ ~ 0.3 GeV obtained from the gap equation. At higher densities widths of e/p remain 
approximately the same. 

To qualitatively characterize this behavior, we evaluate d'^e/dm?, the second derivative 
of e with respect to m at fixed p and m = m*^{p). This quantity describes the stiffness 
of the effective potential along the mass coordinate. It is analogous to the compression 
modulus characterizing the stiffness with respect to the density fluctuation. The results 
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FIG. 11: Second derivative of energy density with respect to nucleon effective mass as a function 
of baryon density. Solid line represents the results of calculation with the parameters from second 
line of Tabled Dashed line corresponds to the same calculation, but with mo = . 

are shown in Fig. ^2 by the solid line. For comparison, the dashed line shows the same 
calculation, but with vanishing bare mass mo . One can clearly see that the stiffness drops 
dramatically at the point of chiral symmetry restoration, where m*^ and ps become small. 
In the case mo = the stiffness vanishes at p = pc and, therefore, the fluctuations diverge. 
Such behavior is expected for a second order phase transition. It is well known that the 
mean-field approximation breaks down in the vicinity of the critical point. However, in a 
realistic case of mo 7^ one can speak only about a significant enhancement of fluctuations 
near the chiral transition point. 

Finally, in Fig. El we present a three-dimensional plot of e/p as a function of m and p, 
calculated for the same set with A = 0.4 GeV. Thick line connecting local minima of this 
surface corresponds to solutions of the gap equation (0). One can easily notice approximate 
symmetry of e/p with respect to reflections m —>■ — m. At low densities one sees two local 



minima 



3l| separated by a barrier (see also Fig. irUj) . At p ^ pc the barrier disappears and 



the surface becomes rather flat along the mass axis. 

It is interesting to note that in the case of exact chiral symmetry (mo = 0) the energy 
functional would be an even function of the effective mass m, i.e. e(— m, p) = e(m,p). 
In particular, at p ^ it has two degenerate minima at m = ±171^. This symmetry is 
simply a particular case of the chiral rotation corresponding to angle n . By using this fact 
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FIG. 12: Energy functional in the coordinates "effective mass-baryon density" calculated within 
the NJL model with A = 0.4 GeV. Thick line shows local minima given by the gap equation. The 
dot marks the global minimum corresponding to the nuclear saturation point. 

one can prove the following theorem, valid also at T 7^ . Namely, the thermodynamic 
potential of any chirally-symmetric theory, Q = —P{'m, T, fi) is an even function of the 
fermion mass m and has an extremum (maximum, minimum or saddle point) at m = . 
Explicit symmetry breaking terms make Q slightly tilted towards m > . All known chiral 
models obey this theorem. In contrast, traditional nuclear models of Walecka type strongly 
violate it. They assign no special significance to the point m = . Such models are tuned 
to the saturation point of nuclear matter and can not be extrapolated far away from this 
point. As an example, in Fig. we show the energy surface predicted by the TMl version 
of the relativistic mean-field model [2^. One can see that the model does not reveal any 
symmetry with respect to m ^ —m . Moreover the minima at m < do not appear at all. 
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FIG. 13: Same as Fig. [El but calculated within the TMl model. 
V. SUMMARY AND DISCUSSION 

In this paper we describe cold nuclear matter using a generalized version of the NJL model 
including additional scalar-vector interaction terms. We formulate the model in terms of 
nucleonic degrees of freedom bearing in mind that normal nuclei are indeed made of nucle- 
ons. We have demonstrated that this model is able to reproduce well observed saturation 
properties of nuclear matter such as equilibrium density, binding energy, compression mod- 
ulus and nucleon effective mass at p = po -The best fit of nuclear properties is achieved with 
the cut-off momentum A ^ 0.4 GeV which is noticeably smaller than usually assumed for 
quark-based models. This is an indication that using nucleonic quasiparticles is only justi- 
fied at low momenta p < A . At higher p the quark degrees of freedom should be included 
explicitly. 

The model predicts two interesting features. First, it reveals a first order phase tran- 
sition of the liquid-gas type occurring at subsaturation densities. This phase transition 
automatically follows from the existence of the nuclear bound state and, therefore, must be 
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present in any realistic model of nuclear matter. Second, the model predicts an approximate 
restoration of chiral symmetry at high baryon densities, p ^ 3po • This feature is also quite 
natural for the NJL model and is expected in any chiral model. But as follows from our 
analysis, this chiral transition does not produce any significant peculiarities in thermody- 
namical quantities. Only in the case of exact chiral symmetry (mo = 0) this phenomenon 
would correspond to a true second order phase transition. It is demonstrated that the chiral 
transition is manifested by increased fluctuations of the order parameter, i.e. the scalar 
condensate. 

We want to stress that our model is able to describe simultaneously the saturation prop- 
erties of nuclear matter and restoration of chiral symmetry at high baryon densities. The 
model exhibits only one first order phase transition of the liquid-gas type at subsaturation 
densities. Restoration of chiral symmetry develops gradually with increasing baryon den- 
sity and does not lead to any phase transition. This confirms conclusions of Ref. where 
we have shown that some popular predictions of the first order chiral phase transition are 
incompatible with phenomenological constraints. 

In conclusion, we have achieved a good description of cold nuclear matter on the basis 
of a relatively simple chiral model. Within this model normal nuclei are interpreted as 
droplets of baryon-rich matter in a phase with spontaneously broken chiral symmetry. The 
gradual approach to a phase with restored chiral symmetry is predicted at baryon densities 
above 3po • This model can be easily extended to nonzero temperatures and finite nuclei. 

Finally we want to emphasize the following. The liquid-gas phase transition discussed 
above can be observed via multiple production of nuclear fragments in nuclear reactions 
at intermediate energies. Nuclear multifragmentation is well established phenomenon and 
many observations reveal trends expected for such a phase transition. The situation with 
the chiral/deconfinement transition in baryon-rich matter is less certain. Despite of many 
interesting suggestions a similar unique phenomenon linked to this transition has not been 
established yet. We believe that effects of the chiral symmetry restoration will be most 
clearly seen in nuclear collisions at energies of a few 10 AGeV, when highest baryon densities 
are expected. But to find them one should be prepared for a thorough study of observables 
sensitive to chiral fiuctuations. 
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